On the extrapolation to ITER of discharges in pressent tokamaks 
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An expression for the extrapolated fusion gain G — Pf us ion/5Phcat (-Pfusion being the total fusion 
power and Phcat the total heating power) of ITER in terms of the confinement improvement factor 
(H) and the normalised beta (/3jv) is derived in this paper. It is shown that an increase in nor- 
malised beta can be expected to have a negative or neutral influence on G depending on the chosen 
confinement scaling law. Figures of merit like H(3n should be used with care, since large values 
of this quantity do not guarantee high values of G, and might not be attainable with the heating 
power installed on ITER. 
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I. INTRODUCTION 

In many tokamak devices discharge scenarios are stud- 
ied with the aim of improving the performance of fu- 
ture reactor experiments over the current design val- 
ues. Essentially two ingredients enter in the optimisation: 
the energy confinement time and the Magneto-Hydro- 
Dynamic (MHD) stability limit, represented by a critical 
pressure. Both energy confinement time and obtainable 
pressure are measured in current experiments and then 
for scaling purposes expressed in dimensionless parame- 
ters. The confinement improvement is given by the so- 
called ii-factor which measures the energy confinement 
time (te) relative to a scaling law (rE,acaling)- The ob- 
tainable volume averaged pressure ((p)) is expressed in 
the normalised quantity /3n 



H = T E /TE, S caling AsT = (A) — 



aB 



(1) 



In the equations above ((3) is the volume average of 
2/iop/B 2 measured in %, a is the minor radius in m, 
B is the magnetic field in T, and I p is the plasma current 
in MA. For the extrapolation one assumes a constant H- 
factor and then uses the confinement scaling to determine 
the confinement time in a next step device. Furthermore, 
it is assumed that the working point of the reactor is at 
the same normalised pressure The H- factors de- 

scribe our imperfect knowledge of the scaling of confine- 
ment, i.e. the confinement in current day experiments in 
some areas of the parameter space scales differently than 
the developed scaling laws suggest. Of course, the use of 
a constant fl-factor is a daring approach to correct for 
this lack of knowledge, but nevertheless can give some 
idea of how confinement could be different in the next 
step experiment. 

Different scenarios for improved performance have 
been proposed (see for an overview These scenar- 
ios include the internal transport barriers (ITB) as well 
as the different scenarios for the improvement of the H- 
mode. The results presented in this paper can in prin- 
ciple be applied to all scenarios, but as examples only 



the H-mode scenarios will be shown The latter scenarios 
are an active area of research with many contributions 
from different machines, for instance, ASDEX Upgrade 
0, DIII-D 0, JET [|, and JT60-U [5]. (Although the 
original reference Q refers to an internal transport bar- 
rier it was shown in a later paper [|| that the ion tem- 
perature profiles follow the same scaling as those of the 
standard H-mode). 

An improvement in confinement or MHD stability 
would allow to operate the next step tokamak experiment 
ITER 0, H| at a higher energy multiplication factor Q, a 
higher fusion power, or a higher bootstrap current frac- 
tion. We define the energy multiplication factor Q and 
the fusion gain G as 



Q 



Pfu 



G 



fusion 



5Phe 



(2) 



where Pf us ion is the total fusion power, P aux is the exter- 
nally applied heating power, and Phcat is the total heating 
power of the plasma, which in steady state is equal to the 
loss power Pi oss - Since one fifth of the fusion power heats 
the plasma, Phcat = Pfuskm/5 + P aU x and there is a direct 
relation between G and Q 



G 



Q 



(3) 



The currently proposed next-step experiment ITER [8| 
is designed to reach Q = 10, and any realistic scenario to 
be tested in this experiment should reach a Q value sig- 
nificantly larger than one. Scenarios at sufficiently large 
Q might be further optimised to reach a higher fusion 
power, or to reach a higher bootstrap current fraction. 
The latter optimisation is known as the hybrid scenario. 
It aims at the extention of the pulse length at similar 
performance as the design value. It is important for the 
results presented in this paper to stress that improve- 
ments in fusion power or bootstrap fraction can only be 
of interest if the the energy multiplication factor is suf- 
ficiently high. In the opinion of the authors a good rep- 
resentation of the extrapolated performance of current 
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discharges towards ITER can only be obtained if one of 
the figures of merit is directly connected with Q. Other 
figures could be used to measure the bootstrap current 
and fusion power. Although Q plays a central role in 
the ITER experiment their is no published simple scal- 
ing that allows a direct assessment of the extrapolated 
value of current discharges. A simple way to obtain a 
rough estimate can be extremely useful in assessing the 
progress made in this large area of research. Finally, we 
would like to stress that our extrapolation formula aims 
at judging the performance in ITER and is not necessar- 
ily applicable for more general purposes. In particular we 
use a fixed size and magnetic field. A reactor design can 
also be optimised through changes in these parameters. 



II. SCALING OF G 



without considering a scaling of this density with plasma 
current) since the density dependence that enters can be 
easily identified through the coefficient a n , which can be 
set to zero to obtain the result for a given fixed density. 
Note that we cannot assume Pi OS s to be constant, since 
discharges at different beta will extrapolate to a differ- 
ent fusion power and, hence, to different plasma heating 
powers. Since for fixed magnetic field and plasma shape 
I p oc qg 5 , with qgs being the safety factor at 95% of the 
plasma radius, and redefining the constant C x to include 
all constant design quantities one obtains 

te = H x r x = H X C X g 9 T— P^ t (8) 

Combining Eq. 1(5)) with Eq. ^ one obtains an expression 
for G. Indicating all quantities of the ITER standard 
scenario with an index S. one can build the ratio 



The fusion power is proportional to 
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where T is the plasma temperature, V is the plasma vol- 
ume, and (<jv) is the cross section for the fusion reactions 
averaged over the velocity distribution. Over the temper- 
ature range of interest (av) oc T 2 , such that the fusion 
power scales with the pressure (p) squared. The power 
loss from the plasma (Pi OS s) is measured by the energy 
confinement time te, and under stationary conditions it 
is balanced by the total heating power 



W 
te 
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2te~ 



heat ) 



(5) 



where W is the stored energy. Combining the Eqs J5J), 
Q, ©, and © one obtains 



G = 



Q 



p 



fusion 



Q 



5P hc 



oc tiTte, 



(6) 



where p — nT was used, and n is the plasma density. 
This is of course the famous tiTte product. 

To proceed we write the confinement time as a product 
of the improvement factor H x over the confinement time 
of an arbitrary scaling law (t x ) 



G _ H x ggss"' 1 Kelt P 
G S H xsq ^P^ p s 



(9) 



Then recalling that P hcat = Pfusion/5G and P fus i on oc 
p 2 oc (/3/v/<fo5) 2 one arrives at 



l+2a p +a/+a„ / a \ l+2a„ 
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Combining the terms containing G we finally derive at 
the desired expression 

1 l+2a p +n;+a„ l+2a p 

G { H x \ f q 95S \ ff3 N ^ 1+a 
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The figures of merit for different scalings can now be de- 
rived directly. Here the explicit expressions are given 
for four different scaling lows. The most commonly 
used IPB98(y,2) indicated by rjf , the L-mode scaling 
ITER89-P denoted by r^, a newly derived scaling from 
Ref. [lj| denoted by Tg, and an electro-static gyro-Bohm 
scaling law derived in [ll| denoted by t§ gb 

rf = 0.I45 J ffH/ p - 93 S 015 ^ a69 n°- 41 M 019 i 1 ' 1 - 39 a - 58 K - 78 

(12) 

r| = 0.048i/ L /;- 85 B - 2 P h -° at 5 n 01 A/ - 5 i? 1 - 2 a - 3 K - 5 (13) 



te = H x t x = H X C X I? 1 B aB P 



loss 



M c 



(7) 

Here C x is a constant, M is the effective mass in AMU, 
and k is the plasma elongation. The exponents a are the 
exponents of the scaling law. For the projection to ITER 
the plasma size (a,R,K) as well as the magnetic field B 
and the effective mass M are assumed to be given by the 
design values. For the density a fixed ratio of the Green- 
wald limit (no) @ will be used, i.e. n oc nQ T I p . In 
our final result it will not be difficult to obtain the result 
at constant density (i.e. the design value of the density 



t% = Om2H c I p 85 B 017 P^ 5 n - 2 Hd 011 R 1 - 21 a - 39 K OS2 
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0.83 D 0.07 D -0.55„0.49^0.14^ 



(14) 



heal 



1.81 0.30 0.75 
(X ri 



(15) 

In the equations above n is the density in units of 10 20 
m~ 3 , R is the major radius in m, a is the minor radius 
in m, M is the averaged ion mass in AMU, K = A/ira 2 
is the plasma elongation, A is the area of the poloidal 
cross section, and -FWt is in MW. The new scaling laws 
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E' r E GB ) have been obtained after designed experi- 
ments have shown a small and possibly absent j3 depen- 
dence of the confinement [TJ El EU ED, EH m contrast 
with the IPB98(y,2) scaling which, when expressed in 
normalised quantities (normalised pressure /3 cx nTB~ 2 , 
normalised Larmor radius p* cx VT/Ba, and normalised 
collisionality v* cx naT~ 2 ) has an unfavourable beta [?J 
dependence 



(16) 



For the electro-static gyro-Bohm scaling, therefore, zero 



beta dependence as well as te cx 
derive 



Bt e cx p-^V 



3 oO, -0.14 



were imposed to 



(17) 



Several papers have pointed out the fact that the absence 
of the beta dependence in te leads at high normalised 
pressure to more optimistic projections for ITER com- 
pared with IPB98(y,2) scaling pjg, El- 

To obtain the scaling expressions for G a standard sce- 
nario must be defined. Here qg§s = 3 and Pns = 1.8 
will be used. The H-factors (Hxs) can be calculated by 
dividing the target confinement time by the confinement 
times of the scaling calculated using Eqs. (fl2"|) . (fl~3|) . (fl"4|) . 
and (fT"5|) . In the latter equations the ITER parameters 
(J p = 15 MA, B = 5.3 T, R = 6.2 m, k = 1.75, n = 10 20 
m" 3 , a = 2 m, M = 2.5, P = 87 MW, t b = 3.68 s) are 
used, yielding Hhs = = 2.2, Hqs = 1-07, and 

Hegb.s — 0.8. One then directly finds 



G — CH x qg 5 Pfj. 



(18) 



with the values of the constant C and the scaling poten- 
tial giving in Table [U From this table it can be seen that 
Pn has a strongly negative effect in the IPB98(y,2) scal- 
ing, and a rather small effect in all the other scalings. It 
is clear from Eq. (fTTj) that G cx (3% occurs only for the 
generic scaling te cx P~ 5 . 

For the derivation of a figure of merit one often ap- 
proximates the coefficients of the scaling law. For a bet- 
ter comparison we can make the similar approximations, 
i.e. a/ = 1, a n = 0, and a p = -2/3 for the IPB98(y,2) 



and a p — —1/2 for all other scaling laws. This yields 



H 3 

G=1Q.8- H 
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7L = 1.24 7c = 5.25 j E gb = 9.375. 



(19) 



(20) 



(21) 



The figure of merit H Pn / q 2 5 can be considered a scaling 
for G if the exponents in the scaling law for the con- 
finement are a>i — 1, a n = 0, and a p — 0, i.e. only 



TABLE I: Values of the constants of Eq. (JTSJ for the different 
scalings. Both coefficients assuming ITER operation at fixed 
absolute density as well as at fixed Greenwald fraction (n = 
0.85ncr) are given. The latter are indicated by the letters 
"Gr". 





C 


X 


Y 


Z 


IPB98(y,2) 


9.62 


3.22 


-1.77 


-1.23 


IPB98(y,2) Gr 


41.15 


3.22 


-3.19 


-1.23 


ITER89-P 


0.892 


2.0 


-1.70 





ITER89-P Gr 


1.11 


2.0 


-1.90 





Cordey 


3.53 


1.82 


-1.72 


0.18 


Cordey Gr 


5.94 


1.82 


-2.20 


0.18 


EGB 


7.41 


2.22 


-1.62 


-0.22 


EGB Gr 


24.52 


2.22 


-2.71 


-0.22 



in the ideal condition if it is assumed that there is no 
power degradation. We note that one can define a fig- 
ure of merit in different ways and that -ff/3/v/<?! 5 can 
be thought of as a combination of good confinement, 
and high fusion power [16| . Also for the generic scal- 
ing (ap = —0.5) the boots trap current in a reactor close 
to ignition scales as H(3n 17| . It is, however, clear that 
H/3jy/q 2 5 does not provide a figure of merrit for G. This 
point is important because it shows that one must be 
careful in using if/Jjv/ggg. From the scaling Eq. (fTTH) 
it follows that the difference with the figure of merit 
iJ/3jv/(?95 i s {^-&Hh/ Pn) 2 ■ The latter quantity is for 
discharges with normal confinement Hh — 1 but high 
normalised pressure /3jy = 3.6 as small as 0.25. This 
makes a large difference in G. The above example with 
Hh = 1, Pn = 3.6, and (795 = 4.2 reaches a value for the 
parameter HhPnI^ that suggest that the ITER tar- 
get (Q = 10) could be reached, whereas in reality such 
a discharges would extrapolate to Q = 1. Not only is 
this value of little interest, it also requires a rather large 
heating power, since the fusion power is four times larger 
and energy multiplication is much smaller. The other 
scaling laws suffer from the same problem, although it is 
less dramatic due to the different exponents of H as well 
as Pn. 



III. DIMENSIONLESS VARIABLES 

Some confusion can arise when considering the scaling 
in terms of the dimensionless parameters P, normalised 
Larmor radius p* = p/a and normalised collisonality v* . 
The scaling of G in terms of dimensionless parameters 
yields [3 



G cx P{Bt e )B cx pi 



1.5+X fll.25+y ,, -0.25+2 p-1-25 



k /3 



v 



(22) 



where X, Y, Z are the coefficient of the p* , j3 and scal- 
ing of confinement as defined in Eq. (|16[) . In the equa- 
tion above the dependence on as well as has been 
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FIG. 1: (color online) Figure of merit HhPn /gls of the 
advanced scenario discharges from ASDEX Upgrade. The 
symbols correspond to different values of q: circles (mag- 
neta) q < 3.5, squares (blue) 3.5 < q < 4.0, stars (black) 
4.0 < q < 4.5, diamonds red 4.5 < q < 5.0, crosses q > 5.0 



even for the electro-static gyro-Bohm scaling. Therefore 
we arrive at the conclusion that for a given design reach- 
ing the beta limit does not help in increasing G. Of course 
the density scaling is more hidden in our approach since 
it is considered to be a design value. One can derive that 
G increases strongly with density for all scalings. 

Finally, it is noted here that for constant heating power 
(3n cx H and all figures of merit have the same form 



Hp. 



N 



4-> 



Av<?9 5 



995 



(24) 



IV. PERFORMANCE DIAGRAM 

Having derived a simple expression that directly al- 
lows to evaluate G, a suggestion is made in this section 
for a diagram that should allow for an easy assessment 
of the extrapolated performance of any discharge. From 
the discussion above it is clear that a better representa- 
tion of the data can be obtained by plotting the scaling 
for G — Q/{Q + 5). against either the scaling of the fu- 
sion power, i.e. (3ff/qg 5 , or the scaling of the bootstrap 
current, i.e. (?95/3at. Since both are of importance it is 
useful to mark the different qg$ values by different sym- 
bols in whatever plot one chooses. This makes that all 
important parameters can be estimated from the same 
graph. 



explicitly added compared with Ref. 11511 . Using the 
scaling G cx /3 1,25+y it was concluded [15| that for the 
IPB98(y,2) scaling (X — —0.9) there is no large bene- 
fit of going to high (3 since Q increases only moderately 
with (3 (/3 35 ), whereas for the energy confinement scal- 
ings that have no beta dependence it would be largely 
advantageous to go to the (3 limit since G cx /3 125 . This 
conclusion seems in disagreement with the results derived 
in this paper, which rather point at a decreasing G with 
(3n for the IPB98(y,2) scaling and a G independent of (3n 
for the other scaling laws. Eq. (|2"2"|) is the correct dimen- 
sionless expression, but it must be noted that the scaling 
with beta G cx /3 125+F holds only at constant normalised 
Larmor radius (p*) as well as collisionality (1/*). The dif- 
ference with the results in this paper is that the results 
are evaluated for a fixed machine size, density, and mag- 
netic field. With these assumptions it is not possible to 
change [3 independently of p* and v*. At fixed density 
the (3 scaling is essentially a temperature scaling, leading 
to changes in the normalised Larmor radius as well as 
the collisionality. Using the scalings of p* and one can 
derive from Eqs. (fT6"|) , (fTT|) , and 



G H ocT 



-1.23 



Gegb oc T 



-0.22 



(23) 



These scalings are consistent with the diagrams of Ref. 
[l5| where G can be seen to decrease with increasing T 
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FIG. 2: Fig. [2] (color online) Scaling for the advanced scenario 
discharges based on the IPB98(y,2) scaling. Symbols reflect 
the 595 values as in figure [T] 
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Figure Q] shows a dataset of advanced scenario dis- 
charges from ASDEX Upgrade in the representation us- 
ing HhPn /<?95 versus ggs/Jjv- In this figure the different 
595 values are indicated with with different symbols (and 
colours in the online version). In the representation us- 
ing H/3n /<795 even the points at highest <?95/3/v reach the 
ITER target. 

Figure [21 shows the scaling derived from Eq. (TT!))) . This 
can be directly compared with Fig.rjjsince the same scal- 
ing law is used. The obtained picture is different in the 
sense that the highest ggs/J/v values no longer reach the 
ITER target for Q. These discharges have only moder- 
ate confinement improvements and high /3jv leading to a 
relatively small (Hh / (3n) 2 ■ In the diagrams Q can ex- 
ceed infinity, which is obviously unphysical. For those 
discharges for which Q > 00, the temperature will rise 
until the fusion cross section no longer scales quadrat- 
ically with T, violating the original assumptions in the 
derivation, and leading to a smaller increase of the fusion 
power with T. This stabilises the solution and leads to 
Q = 00. 

In the diagrams presented so far the external heating 
power is still implicit. A better insight of how much ex- 
ternal heating power is needed to run a certain discharge 
under reactor conditions can be obtained from the dia- 
gram that has the scaling of the fusion power Pf us ion oc 
(/Sjv/ggs) 2 on the x-axis. Because 

-^aux — ^fusion 

IQ one 

obtains 



attainable beta. 



oc 



Q \ Q95 



(25) 



For a fixed auxiliary heating power the relation above de- 
termines a curve in the G = Q/(Q + 5) versus (/Jw/ggs) 2 
diagram. Figure [3] shows the same data as the dia- 
grams before, plotting G — Q/(Q + 5) versus the fusion 
power normalised to the ITER value 2.77 '((3 W / '<Zgs) 2 ■ The 
dashed lines in this diagram are the different values of the 
auxiliary heating power. From left to right I, 2, 4, and 8 
times the ITER design value. This diagram shows that 
high pressure discharges at low Q values would require 
a large amount of installed heating power to be run. It 
is this diagram that carries the largest amount of infor- 
mation and we propose it to be used when representing 
larger datasets of advanced scenario discharges. 

The diagrams presented in this paper are far from per- 
fect, with several effects not properly accounted for: Ra- 
diated power due to Bremsstrahlung as well as dilution 
of the fuel are not properly scaled. The approach with 
constant _ff-factor and (3m is always daring for an extrap- 
olation. If better MHD stability is, for instance, reached 
through current profile shaping, then one should investi- 
gate if such a shaping is extrapolatable to reactor param- 
eters. Also, although (3n is a good scaling quantity for 
ideal MHD instabilities, it does not provide a very good 
scaling for the NTM, which is often found to limit the 
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FIG. 3: (color online) Extrapolated values of Q/(Q + 5) 
as a function of the normalised fusion power 2.77(/3jv/g9s) 2 . 
Symbols reflect the 595 values as in figure [l] The dotted 
lines indicate the amount of external heating (P a ux) necessary 
to run the discharge in ITER. From left to right the curves 
correspond to power levels 1, 2, 4, and 8 times the nominal 
ITER value 



large data sets these diagrams are certainly useful. Also 
they give an idea of what parameters are important when 
developing scenarios. 



V. CONCLUSIONS 

In this paper a scaling for tokamak discharges is de- 
rived that directly measures the fusion gain G = Q/{Q + 
5), and which are consistent with the underlying scaling 
laws. It is shown that (3m does not have a positive in- 
fluence on G, although it does of course extrapolate to a 
larger fusion power. Care is to be taken with figures of 
merit like H(3n /q^- Although this figure of merit does 
measure a combination of good confinement and high fu- 
sion power, the ITER target value of such a quantity does 
not automatically imply a discharge with a sufficiently 
high Q, and might not be attainable with the limited 
heating power that will be installed. A proposal is made 
for a graphical representation in which both the extrap- 
olated Q as well as the fusion power or the bootstrap 
fraction can be directly assessed. 



G 



[1] X. Litaudon, E. Barbato, A. Becoulet, et al., Plasma 

Phys. Contr. Fusion (2004), 46 A19 
[2] O. Gruber, R.C. Wolf, R. Dux, et al., Phys. Rev. Lett. 

(1999), 83 1787 
[3] T.C. Luce, M.R. Wade, P.A. Politzer, et al., Nucl. Fusion 

(2001) , 41 1585 

[4] A.C.C. Sips, E. Joffrin, M. de Baar, et al, Proceedings of 
the 30 th EPS conference on controlled fusion and plasma 
physics, St. Petersbrug Russia 2003, edited by R. Koch 
and S. Lebedev [Europhys. Conf. Abstr. (2003) 27A O- 
1.3A ] 

[5] Y. Kamada, A. Isayama, T. Oikawa, et al, Nucl. Fusion 

(1999), 38 1845 
[6] A.G. Peeters, O. Gruber, S. Giinter, et al, Nucl. Fusion 

(2002) , 42 1376 

[7] ITER Physics expert groups on confinement and trans- 
port and confinement modelling, Nucl. Fusion (1999), 39 
2175. 

[8] M. Shimada, D.J. Campbell, M. Wakatani, et al, IAEA- 
CN-77/ITERP/05 



[9] M. Greenwald, J.L. Terry, S.M. Wolfe, et ah, Nucl. Fusion 
(1988), 28 2199 

[10] J.G. Cordey, K. Thomsen, A. Chudnovskiy, et al, Scaling 
of the energy confinement time with /3 and collsionality 
approaching ITER conditions, To appear in Nucl. Fusion. 

[11] C.C. Petty, J.C. DeBoo, R.J. La Haye et al, Fusion Sci. 
Technol. (2003), 43 1 

[12] JET team, presented by J.G. Cordey, Plasma physics and 
controlled nuclear fusion research, Montreal 1996 (Inter- 
national Atomic Energy Agency, Vienna, 1997), 1, 603. 

[13] C.C. Petty, T.C. Luce, J.C. DeBoo, et al, Nucl. Fusion 
(1998), 38 1183 

[14] D.C. McDonald, J.G. Cordey, C.C. Petty, et al, Plasma 
Phys. Control. Fusion (2004), 46 A215-A225 

[15] C.C. Petty, T.C. Luce, D.C. McDonald, et al, Physics of 
Plasmas (2004), 11 2514 

[16] R.C. Wolf, Plasma Phys. Control. Fusion (2003), 45 Rl 

[17] A.G. Peeters, Plasma Phys. Contr. Fusion 42 B231 
(2000) 



